A STOCHASTIC LAGRANGIAN PROOF OF GLOBAL 
EXISTENCE OF THE NAVIER-STOKES EQUATIONS FOR 
FLOWS WITH SMALL REYNOLDS NUMBER. 



Abstract. We consider the incompressible Navier-Stokes equations with spa- 
tially periodic boundary conditions. If the Reynolds number is small enough 
we provide an elementary short proof of the existence of global in time Holder 
continuous solutions. Our proof is based on the stochastic Lagrangian formu- 
lation of the Navier-Stokes equations, and works in both the two and three 
dimensional situation. 



describe the evolution of the velocity field of an incompressible fluid with kinematic 
viscosity v > 0. One of the (still open) million dollar problems posed by the Clay 
Institute [9] is to show that given a smooth initial data uq the solution to (|1.1|) - (|1.2|) 
in three dimensions remains smooth for all time. 

In two dimensions, the long time existence of (|1.1|) - (|1.2|) is well known 3J. In 
three or higher dimensions, long time existence is known provided a smallness 
condition is imposed on the initial data (see for example [16] for a criterion which 
in some sense is the most general smallness condition). Recently Chemin and 
Gallagher found a (non-linear) criterion on the initial data that guarantees global 
existence of (|l.l|l - (|1.2|) . and does not reduce to a smallness criterion in BMO -1 . 

In this paper we prove global existence of (| 1 . 1 [) — (11 . 2[) provided our initial data 
has small Holder norm. Though global existence under our assumptions can be 
deduced from the Koch-Tataru result, the proof we present here (Section [5]) is 
short, 'elementary' and essentially relies only on the decay of heat flows (Section 
SI), and a stochastic representation of the Navier-Stokes equations using particle 
trajectories (Section [3j see also [5|I14]). 

2. NOTATIONAL CONVENTIONS AND DESCRIPTION OF RESULTS 

In this section we describe the notational convention we use, and state the main 
result we prove. Let L > be a fixed length scale, and X = [0, L}. We define the 
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1. Introduction 



The Navier-Stokes equations 



(1.1) 
(1.2) 



d t u + (u ■ V)u ~ vAu + Vp = 
V • u = 
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Holder norms and semi-norms on X d by 

I, _ L a Hx)-u(y)\ 

x,y£l d F V\ 



u 



u 



C k 



E Lk \ Dm < 



\m\—k 

where D m denotes the derivative with respect to the multi index m. We let C k de- 
note the space of all fc-times continuously differentiable spatially periodic functions 
on I, and C k > a denote the space of all spatially periodic k + a Holder continuous 
functions. The spaces C k and C k,a are endowed with the norms ||-|| cfc and ||-|| fcQ , 
respectively. 

We use I to denote the identity function on R d (or on T d depending on the 
context), and use I to denote the identity matrix. The main theorem we prove in 
this paper is 

Theorem 2.1. Let k ^ 1, a S (0,1) and uq E C k+1,a (X d ) be spatially periodic, 
divergence free and have mean 0. Let R = jj\\uo\\ k+1 a be the Reynolds number of 
the flow. Then 3T = T(k,a,d, j^\\uo\\ k+1 a ) and Rq = Ro(k, a, d) such that for all 
R < Ro the solution u of (jl.l|l - (|1.2|) with viscosity v = j^\\ u o\\k+i a ' *™^ a ^ data 
uq and periodic boundary conditions is in C k+1,a for time T , and satisfies 

(2- 1 ) \\ u T\\ k+ha < lluollfc+La 

We prove Theorem 12. II in Section [5j A few remarks are in order. 



Remark 2.2. Local existence (Theorem IA.5j) combined with the Theorem [2J] im- 
mediately show that for given initial data, we can choose v large enough so that 
(|1.1|) - (|1.2|) have time global C k+1 ' a solutions. Alternately for fixed viscosity, if 
ll M o||fc + iQ, is small enough, Theorems 12.11 and IA.5I again give time global C k+1 ' a 
existence of (fTTT]) - (fTT^|) . 

Remark 2.3. The assumption that uq has mean is not restrictive. First note that 
our boundary conditions imply that J ut is conserved in time. Set u = J uq 
to be the mean velocity. Now if we change to coordinates moving with the mean 
velocity by letting u'(x, t) = u(x + ut, t) — u, then u' solves (|l.lj) - (|1.2[) with mean 
initial data uq — u. Thus the smallness assumption in Remark l2.2l is really smallness 
assumptions on the deviation from the mean velocity. 



Remark 2.4. Theorem 12.11 shows that for some time T, equation (|2.ip holds. Un- 
fortunately our proof does not show that ||itt|| fe+1 a is decreasing in time. 

3. The stochastic Lagrangian formulation 

The Kolmogorov forward equation (or Feynman-Kac formula) [111115] have been 
extensively used to represent solutions of linear parabolic PDE's as the average of 
a stochastic process. In this section we briefly describe here a different approach 
developed in [5ll6]ll4j. which we use to provide a representation of the Navier-Stokes 
equations based on noisy particle paths. 
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Let u : W 1 x [0, oo) — * M. d be some given (time dependent) vector field, and 9 a 
solution to the heat equation 

(3.1) d t + (u • V)0 - vA9 = 

with initial data 9$. We impose either periodic or decay at infinity boundary 
conditions on 9. 

We express 9 as the expected value of a stochastic process as follows: Let W be 
a d dimensional Wiener process, and let X : R d — > IR d be a solution to the SDE 

dX = udt + V2v dW 

with initial data Xo(a) — a. Standard theorjQ [1BJ shows that the flow X is a 
homeomorphism, and as spatially differentiable as u. We let A t denote the spatial 
inverse of the flow map X t . 

Proposition 3.1. If u £ C , 9o £ C 2 then the unique solution 9 of (|3.ip with 
initial data 9 and either periodic or decay at infinity boundary conditions is given 
by 

(3.2) 6 t = m (A t ) 

where E denotes the expected value with respect to the Wiener measure. 

Note that if v = 0, then Proposition 13. II is nothing but the method of charac- 
teristics. If v > 0, this can be interpreted as solving along random characteristics, 
and then averaging. Notice also that the Wiener process s/2vWt is the natural one 
to consider here, as it's generator is v/\. 

The reason we use the representation (|3.2p and not the Kolmogorov forward 
equation is because the Kolmogorov forward equation in it's natural setting involves 
final conditions, and not initial conditions. Thus the standard method employed 
by probabilists is to make a t = T — s substitution [TU]. The process obtained 
in this manner will have the same one dimensional distribution as the process 
At above, however spatial covariances and gradients of the two processes will in 
general be different. The stochastic representation of the Navier-Stokes equations 
we describe below involves spatial gradients of the flow map A, and for this reason 
our representation will not be valid if we use the Feynman Kac formula. 

We now use Proposition 13 . 1 1 to represent the solution to the Navier-Stokes equa- 
tions as the expected value of a system that is nonlinear in the sense of McKean. 
The essential idea is to find a representation of the Euler equations involving parti- 
cle trajectories jj, and then add noise and average as in Proposition l3.1l (as opposed 
to attempting to use the Kolmogorov forward equation). 

Theorem 3.2. Let v > 0, W be an n- dimensional Wiener process, k 1 and 
uq £ C k+1,a be a given deterministic divergence free vector field. Let the pair u, X 
satisfy the stochastic system 

(3.3) dX t =u t dt+ V2u dW t 

(3.4) A t = X^ 1 

(3.5) u t = EP[(V* At) (u o A t )\ 



See also [5ll6lll4| for an elementary proof for flows of the type we consider here 
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with initial data 

(3.6) X(a,0) = a. 

We impose boundary conditions by requiring u and X — I are either spatially peri- 
odic, or decay at infinity. Then u satisfies the incompressible Navier-Stokes equa- 
tions (|l.ip ~ (|1.2p with initial data Uq. 

Here P in equation (|3.5p denotes the Leray-Hodge projection onto divergence 
free vector fields [2] ■ We remark that (|3.5p is algebraically equivalent to 

(3.7) uh = E[(VI t )«„] o At 

(3.8) U( = -A" 1 Vxki 



and (|3.5p can be replaced with (|3.7p - (|3.8p in Theorem l3.2l Note that (|3.8p is exactly 
the Biot-Savart law. When v = 0, equation (|3.7p reduces to the well known vorticity 
transport for the Euler equations [2J, and in this case (|3.3p - (l3.6p (or equivalently 
the system (|3.3p . (|3.4p . (|3.6|) — (|3.8|) ) are exactly a Lagrangian formulation of the 
Euler equations [4]. 

We do not prove Proposition ^ . 1 1 or Theorem 13 . 21 here, and we refer the reader to 
[5]JT4"] instead. For a generalization of Proposition 13 . II where the diffusion matrix is 
not spatially constant we refer the reader to [6].fT4] . 



4. Decay of heat flows 

In this section we prove a decay estimate for solutions to the heat equation 
with an incompressible drift. Our first estimate is an L°° — > L°° estimate that 
is independent of the drift. A more general L 1 — > L°° version of this estimate 
appeared for example in [7] and [8] . We provide a proof that follows the proof in 
[7] and keeps track of the dependence of the constants on viscosity and our length 
scale L. 

Lemma 4.1. Let u G C 1 ([0, T},T d ) be divergence free, and 9 be a solution to the 
equation (|3.ip with initial data 9q. If 9q is spatially periodic, mean 0, and the 
dimension d~^Z, then there exists an constant c = c(d) such that 

PtWco < 7^72 Woo 

Proof. Let <p be mean zero and periodic, p ^ anc j c _ c(d,p) be a constant that 
changes from line to line. Then the Holder, Poincare and Sobolev in inequalities 
give 

V = yv / v 2p ~ i)/p 

^ WPWl 1 IIVIIl(2p-i)/(p-d 

<cL (4p-d-2)/2 P ||^||l/p|| V ^||(2p-l)/ P 

If we set q — 2^~i this gives 

||V^||i 2 ^cL^ 4 )/ 2 ||^|| 2 L t 9 ||^||Z? 
Now let 9' and 8" to be solutions of (13.111 with initial conditions 8„ and 6i 



J ami w Q 

respectively. Integrating p.ip immediately shows that / 9' and J 9" are conserved. 
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Since 9' and 8" are of constant sign, this means that and are conserved 

in time. Finally, the maximum principle implies that 9' 9 9" , and hence H^H^i 
is nondecreasing in time. 

Thus multiplying (13. 1|) by 9 and integrating over I d gives 

d t \\9\\% = -2v\\V9\\ 2 L2 ^ -cvL^- 4 y 2 \\9\\- L ?\\9\\ 2 L + 2 q ^ -cvL^ d -V/ 2 \\9 Q \\l?\\9\\ 2 + q '■ 

Dividing by and integrating in time gives 

L 2/q-d/2 



Let Vt(u) be the solution operator of (|3.1| . The above estimate shows 

L 2/q-d/2 



(^) 1/<? 



Since u is divergence free the dual operator V*(u) — V(—u), and hence satisfies the 
same bound. Thus 







L^L 2 


\n\ 


L 2 ^L~ 






\l 1 ->l^\ 


\v;\ 


l(L~)-^L2 




\\Vt\ 




\n\ 





Hence 



Ti/q-d Ti/q 
ll llL»^C—- 7 -||0 o || il ^ C - 



•(^)2/gll "Nil ^ (vt) 2 /l* 

Finally, p ^ is the same as q ^ | , and choosing q = 4 concludes the proof. □ 

Remark 4.2. When d = 2, p ^ ^2 needs to be replaced with p > ^4^, and hence 
the above proof will show that for any e > 0, 

Halloo < (^)(d+e)/2 Halloo 

The formulation of the Navier-Stokes equations (|3.3|) - (|3.6j) involves recovering 
the velocity from the inverse flow map via a singular integral operator (either Biot- 
Savart or Leray-Hodge) . The unboundedness of singular integrals on L°° (see [T5] ) 
causes Lemma |4. II to be insufficient for our purposes. We now extend Lemma l4.1l to 
Holder spaces for use in our global existence proof. Using the stochastic flows from 
[5j.fT4] we obtain a Holder estimate for solutions of (|3.1|) in an elementary manner. 

The usual PDE methods [17] provide Holder estimates that grow exponentially 
in time. The estimate we provide here will in general grow exponentially in time, 
however decays in time when the viscosity is large, or drift U is small. 

Lemma 4.3. Let d ^ 3 and u E C fc+1,a ([0, T],l d ) be divergence free and define U 
by 

(4.1) U= sup |K|| fe+1 

te[o,T] 
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Let 80 € C(T d ) have mean 0, and 8 satisfy equation (|3.1|) with initial data 80 . Then 
j-, d, k, a) and a constant c — c{^- 

( L d /Ut\ a ^ 



there exists T' = T'(j-, d, k, a) and a constant c = c(^-, d, k, a) such that 



\W\\k+l,a ^ c y^ t y/2 t [j; J J ll"Ullfc+l,a 
holds for all t € [0, T'] . If d = 2, the above estimate is still true if we replace ^ d/2 

rd+e 

with {l/t e ) ( d+e )/2 for any e > 0. 

Remark 4.4. Note that the growth term is independent of the viscosity, and the 
decay term is independent of the drift u. 

Proof. We present the proof for d ^ 3. The d = 2 case will then follow by replacing 
d with d + e. Define X, A by equations ([3.3p and (|3.4[) respectively. From [511T4] 
and uniqueness of strong solutions to (j3.ip we know 

8 t = E0„ ° A t . 

Let I — A — I be the Lagrangian displacement. First notice that if / S C Q then 
Lemma IA.4I shows 

(4-2) \f°M a ^c\f\ a (j- 

Now, let m a multi index with 1 ^ \m\ ^ k. We note that D m 8 t is a sum of terms 
of the form 



(4.3) D n o | A II DTHit and 



A 4 



l<i< \n\ 



where n^'s are multi indices with \rii\ ^ 1 and |n| + |nj| = \m\. By Proposition 
13.11 we know that E*D n 0o\ A satisfies (|3.1[) with initial data D n &o, and hence by 
Lemma |4. II we know 

\\E[D-8 ]oA t \\ Lx ^ rL " 



{vt) d l 2 ' 

Thus using Lemma [A.41 inequality l|4.2p we have 

(4.4) ll E ^o]oA t || a <c(^+(^ 

Using (|4.4p and Lemma [A.41 we bound the remaining terms of (|4.3p . concluding 
the proof. □ 

5. Global existence 

In this section we prove Theorem 12. II We start with a Lemma involving bounds 
for the Leray-Hodge projection. 

Lemma 5.1. Let k ^ 1, and A,v E (jk+x,a ^ e sttc /j \JA, v are spatially 
periodic. There exists a constant c — c(d, a) such that 

M(y*A)v}\\ k+1 , a <c\\VA\\ k Jv\\ k+ha 



GLOBAL EXISTENCE FOR NAVIER-STOKES WITH SMALL REYNOLDS NUMBER. 7 



Proof. Since P vanishes on gradients, we can 'integrate by parts' to avoid the loss 
of derivatives. Note 

P[(V*tt)«] = P[V(it ■ v) - (y*v)u\ = -P[{V*v)u] 

for any u, v € C . Thus we have 

0iP[(V*A)«] = P[(V*A)c^] - P[(V*«)cU] 

Since P is Calderon-Zygmund singular integral operator, it is bounded on Holder 
spaces [TIIT5] . Finally note that the right hand side only depends on first derivatives 
of A and v, and the lemma follows by taking Holder norms. □ 

We now prove Theorem 12. 11 We restate it here for the readers convenience. 

Theorem (j2~Tj) . Let k > 1, a € (0, 1) and u € C k+1 ' a (l d ) be spatially periodic, 
divergence free and have mean 0. Let R = |r||wo|| fc+1 a be the Reynolds number of 
the flow. Then 3T = T(k,a,d, -^\\uo\\ k+1 a ) and Ro = Ro(k,a,d) such that for all 
R < Rq the solution u of (jl.l|l - (|1.2|) with viscosity v — 4||uo|| fe , i a , initial data 
uq and periodic boundary conditions is in C k+1 ' a for time T, and satisfies 

(EH) Kll*+l,a < H«0lUl,a 

Proof. We assume that d ^ 3. The d — 2 case follows similarly by replacing 
d with d + e. Let C, 8 be the dimensional constants in Theorem IA.5I We let 
U = C\\uo\\ k+1 a , and choose T such that 2£ < S. By Theorem I A. 5 [ there exist a 
pair of C k+1 ' a functions X, u : [0, T] — > L d which are the unique (strong) solution 
to l(3T3 |) -([3l] l . Recall that ||u t || fc+1)Q < U for all t e [0,T]. 
From equation (|3.5p we see 

ut = EPw o A t + EP (V*t) u oA t . 

Let c = c(k, a, d) be a constant that changes from line to line. Applying Lemma 
15.11 and Lemma [A. 41 to the second term we have 

ht\\ k+ha < c\\Eu o A t || fc+1;a + cE||V*^|| fcia ||uo o A t \\ k+l a 
and hence by Lemma l4.3l we have 

J ll U o||fe + i, Q - 

Minimizing ^^ d/2 + (■^■) Q in time shows that the minimum value is attained at 
t = £-kj{ d /( 2a + d ) i an d the minimum value is c R ad /( 2a + d ) . Thus we can choose i? 
small enough to ensure t < T and equation (|2.ip is satisfied. □ 

Appendix A. Bounds for the Lagrangian displacement 

In this section, we prove bounds on ||VX — I|| fc Q . The estimates proved here 
are elementary, and are taken directly from [131114] , We reproduce them here for 
completeness and the readers convenience. 

We remark that the estimates provided here were used in jT3l[14] to prove local 
existence for the system (|3.3j) - (|3.6j) . As the local existence proof is a little lengthier, 
we do not reproduce it here. 



( L d fUt\ 
k+1 > a ^ C {(vt) d /z + {l) 
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Lemma A.l. Let X be a Banach algebra. If x E X is such that \\x\\ ^ p < 1 then 
1 + x is invertible and || (1 + a;) -1 1| ^ T^p' Further if in addition \\y\\ ^ p then 

\\(l + X )-l -(l + y)-l\\<^-j-^\\ X -y\\ 

Proof. The first part of the Lemma follows immediately from the identity (1 + 
x)^ 1 = Yl(~ x ) n - The second part follows from the first part and the identity 

(1 + x)- 1 - (1 + y)- 1 = (1 + xy^y -x)(l + y)- 1 . □ 
Lemma A. 2. Ifk^l, then there exists a constant c = c(k,a) such that 

\\f°9\\ k , a <4f\\ kta (l + ||V<7|| fe _ 1)a ) fe+a 

and 

11/ o gx - f o g 2 \\ Ka < c\\Vf\\ kia (l + ||V5i|| fc _ lia + ||V 92 || fe _ 1)a ) &+1 • 

• llfli - flail fc, Q - 

The proof of Lemma lA. 21 is elementary and not presented here. 
Lemma A.3. Let X t ,X 2 E C k+1 ' a be such that 

ll V ^l -%,a < d < 1 aTld ll V ^2-I||fe !Q <d<l. 

Lei Ai and A2 oe i/ie inverse of X\ and X2 respectively. Then there exists a constant 
c = c(k, a, d) such that 

\\Ai-A 2 \\ kta ^c\\X l -X 2 \\ kia 

Proof. Let c = c(k, a,d) be a constant that changes from line to line (we use this 
convention implicitly throughout this paper). Note first VA — (VX) _1 o A, and 
hence by Lemma I A. II 

||VA|| C „ < \\{VX)- l \\ co <c. 
Now using Lemma [A. II to bound ||(VX) _1 || Q we have 

||VA|| Q = H(VX)- 1 o A\\ a < IKVX)- 1 ^ (1 + ||VA|| C0 ) < c 

When k ^ 1, we again bound ||(VX) _1 || fe a by Lemma |A~T1 Taking the C k ' a norm 
of (yxy 1 o A we have 

\\VA\\ ka ^\\(VX)-\ a (l + ||VA|| fc _ l!a ) fe . 

So by induction we can bound || VA|| fc by a constant c = c(fc, a, d). The Lemma 
now follows by applying Lemma lA. 21 to the identity 

A x - A 2 = {Ax 0X2-1)0 A 2 

= (Ax oX 2 - Ax 0X1)0 A 2 . □ 

Lemma A.4. Let u E C([0, T],C k+1,a ) and X satisfy the SDE (US) with initial 
data ()3.6p . Lei A = X — I and U = sup t ||u(i)||j. +1 a - Then there exists T = 
T(k, a, j^) and c = c(k, a, ^r-) such that for t ^ T 

cTJt cllt 
HVA(t)|| fc)Q < — and ||V*(t)|| fc)Q *S — . 

hold almost surely. 
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Proof. From equation (|3.3j) we have 

X(x,t)=x + [ u(X(x,s),s)ds + V2vW t 
Jo 

(A.l) =^ WX{t) = I + I (Vu)oX-VX. 

•/ o 

Taking the C° norm of equation (|A.1|) and using Gronwall's Lemma we have 

||VA(*)|| C o - ||VJC(i) - 7|| c0 s$ e ut ' L - 1. 
Now taking the C < a norm in equation (|A.1[) we have 

l|VA(t)|| fc , a < cj* ||Vt*|| fc)Q (l + ||VA|| fe _ liQ ) fc+ " (l + ||VA|| fcia ) . 

The bound for ||VA|L a now follows from the previous two inequalities, induction 
and Gronwall's Lemma. The bound for a then follows from Lemma [A. 31 

We draw attention to the fact that the above argument can only bound VA, and 
not A. Fortunately, our results only rely on a bound of VA. □ 

We conclude this appendix by stating a slightly modified version theorem which 
appeared in [13] . The only modification we make is that we trace the dependence 
of the constants in [T5] to dimension less quantities, instead of absolute ones. The 
proof that appeared in [13] goes through verbatim. 

Theorem A. 5. Let k J? 1 and uq £ Qk+i,a y e divergence free. There exists 
absolute constants S = S(k,a,d) and C = C(k,a,d) such that for U — C||uo|| fe+1 a , 
and any T such that ^f- < S there exist a pair of functions a pair of functions 
A, u G C([0, T],C k+1 ' a ) such that u and X = I + A satisfy the system (|3~3"j) - (|3TB|) . 
Further for all t £ [0, T] we have ||Mt||fe +1 a ^ U . 
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